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Abstract 

A design study is currently in progress for a third-generation gravitational-wave (GW) detector called 
Einstein Telescope (ET). An important kind of source for ET will be the inspiral and merger of binary 
neutron stars (BNS) up to z ~ 2. If BNS mergers are the progenitors of short-hard 7-ray bursts, then some 
fraction of them will be seen both electromagnetically and through GW, so that the luminosity distance and 
the redshift of the source can be determined separately. An important property of these 'standard sirens' is 
that they are self-calibrating: the luminosity distance can be inferred directly from the GW signal, with no 
need for a cosmic distance ladder. Thus, standard sirens will provide a powerful independent check of the 
ACDM model. In previous work, estimates were made of how well ET would be able to measure a subset of 
the cosmological parameters (such as the dark energy parameter wq) it will have access to, assuming that 
the others had been determined to great accuracy by alternative means. Here we perform a more careful 
analysis by explicitly using the potential Planck cosmic microwave background data as prior information for 
these other parameters. We find that ET will be able to constrain wq and w a with accuracies Au>o = 0.099 
and Aw a — 0.302, respectively. These results are compared with projected accuracies for the JDEM Baryon 
Acoustic Oscillations project and the SNAP Type la supernovae observations. 

PACS numbers: 98.70.Vc, 98.80.Cq, 04.30.-w 
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I. INTRODUCTION 

In the past decade, various observations, including type-la supernovae (SNIa) the temper- 
ature and polarization anisotropics power spectrum of the cosmic microwave background (CMB) 
radiation [2|, the baryon acoustic oscillations (BAO) peak in the distribution of Sloan Digital Sky 
Survey luminous red galaxies 3], and weak gravitational lensing [4], have all suggested that the 
present Universe is undergoing an accelerated expansion. A possible explanation would be the 
presence of a fluid called dark energy, which should have positive density but negative pressure (for 
a review, see Understanding the physical character of dark energy, assuming it exists, is one 
of the main challenges of modern cosmology. A key question is then how well we will be able to 
differentiate between various dark energy models by measuring the dark energy equation of state 
(EOS) and its time evolution. 

Currently, among the main methods to determine the dark energy EOS are observations of SNIa, 
CMB, and large-scale structure. The capabilities of these methods will be improved significantly 
in the near future However, we note that all these methods are based on the observations 

of various electromagnetic waves. In addition to these electromagnetic methods, the observation 
of gravitational waves (GW) will provide a new technique, where gravitational-wave sources, in 
particular inspiraling and merging compact binaries, can be considered as standard candles, or 
standard sirens g]. In the case of ground-based detectors, the idea is to use binaries composed of 
two neutron stars (BNS), or a neutron star and a black hole (NSBH). These are hypothesized to 
be at the origin of short-hard 7-ray bursts (shGRBs). In many cases it is possible to identify the 
host galaxy of a shGRB and determine its redshift. From the gravitational- wave signal itself one 
would be able to measure the luminosity distance in an absolute way, without having to rely on a 
cosmic distance ladder: standard sirens are self- calibrating. 

The use of GW standard sirens to measure the Hubble constant with a network of advanced 
ground-based detectors has been studied by Nissanke et al. jjj, and with LISA (using extreme mass 
ratio inspirals) by MacLeod and Hogan [101 ] . Supermassive binary black holes may be useful to 



study dark energy with L 
cosmology and gravity [17 



SA 
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161 ]; more generally, they can constrain alternative theories of 



20|. Observations of BNS events with the Big Bang Observer would 



also allow for dark energy studies 



21|. 



Currently a third-generation ground-based observatory called Einstein Telescope (ET) is un- 
dergoing a design study [22]. The latter would be able to see BNS inspirals up to redshifts of 
z ~ 2 and NSBH events up to z ~ 8, corresponding to millions of sources over the course of several 
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years, some fraction of which will have a detectable electromagnetic counterpart (e.g., a shGRB). 
Sathyaprakash et al. have investigated how well cosmological parameters could be determined with 
ET assuming 1000 'useful' sources 



23j. Among the parameters which ET will have access to are 

(H ,tt m ,tt k , w ,w a ), (1) 

where Hq is the Hubble parameter at the current epoch, Q, m the density of matter by the critical 
density, Q,k a parameter related to spatial curvature, and wq and w a parameters determining the 
dark energy EOS and its time evolution (see below for more precise definitions). ET will not be able 



to arrive at a completely independent measurement of all these parameters at once. In [23J it was 
assumed that, e.g., all parameters except wq had been measured by other means (electromagnetic 
or GW) and could be assumed known with arbitrary accuracy for all practical purposes. Here we 
continue this study in more depth, with a focus on the dark energy parameters wq and w a . Instead 
of assuming the other parameters to be exactly known, we will use the predicted CMB prior from 
Planck. CMB measurements give accurate values for Hq, £l m , fifc, but have large uncertainties 
in wq and w a . Heuristically, imposing this prior effectively 'fixes' the values of Hq, Q m , Q^- To 
measure wq, w a with GW standard sirens is then an important check of the values obtained through 
electromagnetic means. 

The outline of the paper is as follows. In Sec. HU we recall the basics of using short-hard 7-ray 
bursts as standard sirens in potential ET-GW observations. We then discuss the determination of 
the dark energy parameters by the ET-GW method alone, after which we impose the Planck CMB 
prior. In Sec. IIIIl we discuss the capabilities of the JDEM BAO project and of the SNAP SNIa 
project, and a comparison as well as the potential combination with the ET-GW method is given. 
In Sec. IIVI we conclude with a summary of our main results. 



II. SHORT-HARD 7-RAY BURSTS AS A KIND OF STANDARD SIRENS 
A. The expanding Universe and the dark energy 

We will work with the Friedmann-Lemaitre-Robertson- Walker universes, which are described 

by: 

ds 2 = _ dt 2 + a 2 {t) I Jf_ + r 2 d Q2 + r 2 gin 2 ^2 j ( (2) 

where t is the cosmic time, and (r, 6, (f>) are the comoving spatial coordinates. The parameter 
k = 0,1,-1 describes the flat, closed and open universe, respectively. The evolution of the scale 
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factor a(t) depends on the matter and energy contents of the Universe. Within general relativity, 
the equations governing this evolution are 

a\ 2 _ „ 2 SnGptot k d 4vrG 

- = H = 2 ' _ = ^-(Ptot + 3ptot), (3) 

a J 3 gt a 3 

where ^j t and are the total energy density and pressure in the Universe, and H is the Hubble 
parameter. Since in this paper we are mainly interested in the later stages of the evolution of the 
Universe, where the radiation component can be ignored, we only take into account baryonic matter, 
dark matter, and dark energy. The baryon and dark matter are both modeled as pressureless dust. 
We will assume that the EOS of a dark energy component is responsible for the recent expansion 
of the Universe, which should be determined by observations. In this paper, we shall adopt a 
phenomenological form for the equation-of-state parameter w as a function of redshift z: 



(z) = Pde/Pde = W + W a (l - a) + O [(1 - a) 2 ] 

~ W + Wa— — ■ (4) 

L ~\~ Z 



This is the so-called Chevallier-Polarski-Linder form [241] .which has been adopted by many authors, 
including the DETF (Dark Energy Task Force) group In the present epoch where z ~ 0, we 
have w ~ wq. w a describes the evolution of w to next-to-leading order in (1 — a). Since we are 
mostly interested in the later stages of the Universe's evolution, higher order terms will be ignored. 

The evolution of dark energy is determined by the energy conservation equation pd e + 2>H{pd e + 
Pde) = 0. By using the EOS of dark energy, Eq. we find that 

Pde = Pdefl x E{z), (5) 
where pd e ,o is the value of p^e at z = 0, and 

E( Z ) = (1 + z f(l+™0+W a ) e -3WaZ/(l+z)_ ( 6 ) 

Using Eq. (|3|), the Hubble parameter H then becomes 

H(z) = H [n m (l + zf + O fe (l + zf + (1 - n m - n k )E(z)} 1/2 . (7) 

In this expression, Q m = 8TrGp rrij o/3HQ is the density of matter (baryon as well as dark matter) 
relative to the critical density, and £1^ = —k/H^ is the contribution of the spatial curvature. Ho 
is the Hubble parameter at the present epoch. Throughout this paper, we shall adopt a fiducial 
cosmological model with the following values for the parameters |2|]: 

w = -1, w a = 0, £l b hl = 0.02267, « c /ig = 0.1131, n k = 0, h = 0.705, (8) 



where ho = -ffo/(100kms _1 Mpc -1 ). The other parameters are obtained as Q m = Qb+£l c = 0.2736, 
f^e = 1 — Q m — fifc = 0.7264. In Sec. Ill Dl the CMB prior for the dark energy determination will 
be discussed, where the perturbation parameters A s and n s (the amplitude and spectral index of 
primordial density perturbations, respectively) and the reionization parameter r (the optical depth 
of reionization) are also needed. In our fiducial model, we take these to be [2j 



A s = 2.445 x 1(T 9 , 



0.96, t = 0.084. 



(9) 



To conclude this subsection, we state the expression for the luminosity distance d>L of the 
astrophysical sources as a function of redshift z (see, e.g., 



d L (z) = (! + *)< 



|*| 

z dz' 



1/2 ■ 

' sm 



\ k \ 1/2 io ( ^ <0) ' 

(n k = o), 



It 



H(z') 

sinh 



1/2 pz dz' 
Hjz 7 ) 



l*r/c 



(10) 



(n k > o), 



where |/c| 1//2 = Hoy/\£lk\- This formula will be used frequently in the subsequent discussion. 



B. Short-hard 7-ray bursts and gravitational waves 

Current observational studies of dark energy strongly rely on standard candles, i.e., sources for 
which the intrinsic luminosity is assumed to be known within a certain tolerance so that they can 
be used to determine luminosity distance. A widely used standard candle is the type-la supernova 
(SNIa) . The intrinsic luminosity of distant SNIa needs to be calibrated by comparison with 

different kinds of closer-by sources, leading to a 'cosmic distance ladder'. This is not the case with 
GW standard sirens. As pointed out by Schutz, the chirping GW signals from inspiraling compact 
binary stars (neutron stars and black holes) can provide an absolute measure of distance, with 
no dependence on other sources [8]. The GW amplitude depends on the so-called chirp mass (a 
certain combination of the component masses) and the luminosity distance. However, the chirp 
mass can already be measured from the signal's phasing, so that the luminosity distance can be 
extracted from the amplitude. 

Before discussing standard sirens in more detail, let us first recall some basic facts about the 
gravitational radiation emitted by inspiraling compact binaries. Gravitational waves are described 
by a second rank tensor h a p, which, in the so-called transverse-traceless gauge, has only two 
independent components h + and h x , h xx = —h yy = h+, h xy = h yx = h x , all other components 
being zero. A detector measures only a certain linear combination of the two components, called 
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the response h(t), which is given by (see, e.g., [261 ]) 



h(t)=F+(e,^)h+(t) + F x (9,<l>,il))h x (t) i (11) 

where F + and F x are the detector antenna pattern functions, ip is the polarization angle, and (9, cj>) 
are angles describing the location of the source on the sky, relative to the detector. In general these 
angles are time-dependent. In the case of Einstein Telescope, binary neutron star signals can be in 
band for hours, but almost all of the signal-to-noise ratio will be accumulated in the final minutes 
of the inspiral process. With LISA, Doppler modulation due to the orbital motion, as well as 
spin precession, will allow for accurate determination of the angular parameters (see, e.g., [27( and 
references therein), but this is unlikely to happen for BNS (or NSBH) signals in ET with Doppler 
modulation due to the Earth's rotation. Nevertheless, some improvement in parameter estimation 
can be expected, which for simplicity we do not take into account here. In the sequel, (9, <f>, ip) will 
be considered constant. 

Consider a coalescing binary at a luminosity distance di, with component masses mi and m,2- 
Write M = mi + m% for the total mass and 7] = mim^/M 2 for the symmetric mass ratio, and 
define the 'chirp mass' as M c = Mr] 3 / 5 . For sources at cosmological distances, what enters the 
waveform is the observed chirp mass, which differs from the physical chirp mass by a factor (1 + z): 
-Mc,obs = (l + a;) -M CiP hy S - Below, M. c will always refer to the observed quantity. To leading order 
in amplitude, the GW polarizations are 

h+(t) = 2M 5 J 3 d L 1 (l + cos 2 (l))uj 2/3 (t -t) cos[2$(t - t;M, 77) + $ ], (12) 
h x (t) = AMl^dl 1 cos(i)uj 2/3 (t - t) sin[2$(t - t; M, rj) + $ ], (13) 

where l is the angle of inclination of the binary's orbital angular momentum with the line-of-sight, 
oj{pQ — t) the angular velocity of the equivalent one-body system around the binary's center-of- 
mass, and <3?(to — t', M, rj) the corresponding orbital phase. The parameters to an d are constants 
giving the epoch of merger and the orbital phase of the binary at that epoch, respectively. The 
phase $ has been computed perturbatively in the so-called post-Newtonian formalism (see [281 ] and 
references therein). Since we will mostly be concerned with binary neutron stars, spi n will not be 
important, in which case the phase is known up to 3.5PN in the usual notation [29|, and this is 
what we will use here. 

During the inspiral, the change in orbital frequency over a single period is negligible, and it is 
possible to apply a stationary phase approximation to compute the Fourier transform H(f) of the 
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time domain waveform h(t). One has 

H{f) = Ar V& ^V [i(2ir/to - vr/4 + 2^(//2) - <^ (2>0) )] , 
where the Fourier amplitude A is given by 

A = J~ L \j F U l + cos 2 (0) 2 + F 2 Acos 2 (l)^tt~ 7 / 6 M 5 J 6 . 



(14) 



(15) 



The functions ^ and </>(2,0) are 



</>(/) 



^(2,0) 



"^0 + 



256?y 



tan 



i=0 

2 cos(i)i ? x 



/3 



(16) 
(17) 



(l + cos 2 (0)i^ 

The parameters tpi can be found in [26|]. %(/) is taken to be zero outside a certain frequency 
range. The upper cutoff frequency is dictated by the last stable orbit (LSO), which marks the end 
of the inspiral regime and the onset of the final merger. We will assume that this occurs when the 
radiation frequency reaches / upP er = ^Ibso-, with fiso = l/(6 3//2 27rM b s ) the orbital frequency at 
LSO, and M Q b s = (1 + z) M p ^ ys the observed total mass. 

In this paper we shall focus on the observation of GW sources by the Einstein Telescope, a 
third-generation ground-based gravitational-wave detector. Although the basic design of ET is 
still under discussion, one possibility is to have three interferometers with 60° opening angles and 
10km arm lengths, arranged in an equilateral triangle 
functions are: 



22( | . The corresponding antenna pattern 





6,V) 


V3 
2 


-(1 + cos 2 (<9)) cos(20) cos(2-i/0 - 


cos(0) sin(2</>) sin(2^) 






V3 
2 


-(1 + cos 2 (6»)) cos(2^) sin(2^) + 


cos(#) sin(2</>) cos(2?/>) 








(0,^ + 2vr/3,^), 




F?Ue,< 


6,V) 




(0,0 + 4vr/3,^). 





(18) 

The performance of a GW detector is characterized by the one-side noise power spectral density 
Sh(f) (PSD), which plays an important role in the signal analysis. We take the noise PSD of ET 
to be 

1 + b\x + b2X 2 + b 3 x 3 + 64X 4 + 65X 5 + bQX 6 



30] [45] 

S h (f) 



So 



x Pl + a x x iJ1 + a 2 



P2 



1 + C\X + C2X 2 + C3X 3 + C4X 4 



(19) 
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where x = f/fo with /o = 200Hz, and So = 1-449 x 10 52 Hz 1 . The other parameters are as 
follows: 

Pl = -4.05, p 2 = -0.69, 
a x = 185.62, a 2 = 232.56, 
bt = 31.18, b 2 = -64.72, b 3 = 52.24, b A = -42.16, 6 5 = 10.17, 6 6 = 11.53 

ci = 13.58, c 2 = -36.46, c 3 = 18.56, c 4 = 27.43. (20) 

For data analysis proposes, the noise PSD is assumed to be essentially infinite below a certain 



lower cutoff frequency /i OW er (see the review [26f]). For ET we take this to be /i OW er = 1 Hz. 

The waveforms in Eq. (114p depend on the seven free parameters 
(In .M c , In//, to, cos(t), ip, In d^); note that for 'useful' events the sky position will be known. In 
order to deal with the parameter estimation, throughout this paper, we employ the Fisher matrix 
approach 



311 ] . Comparing with the Markov chain Monte Carlo (MCMC) analysis, the Fisher 
information matrix analysis is simple and accurate enough to estimate the detection abilities of 
the future experiments. In the case of a single interferometer A (A = 1,2,3), the Fisher matrix is 
given by 

A^ = (n^,-Hf\ nl A) = dn^(f)/d Pl , (21) 

where is the output of interferometer A, and the pi denote the free parameters to be estimated, 
which are 

(In M c , In T], t , $ , cos(t), -0, In d L ). (22) 

The angular brackets denote the scalar product, which, for any two functions a(t) and b(t) is 
defined as 

(l>i)B< ^ W)^-(f)itf) l m 

J /lower 1 b h{J) 

where a and b are the Fourier transforms of the functions a(t) and b(t). The Fisher matrix for the 
combination of the three independent interferometers is then 

3 



A, ; (24) 



A=l 



The inner product also allows us to write the signal-to-noise ratios A = 1, 2, 3 in a compact 
way: 

p(A) = ^ H (A) jH (A)^ ( 25) 



The combined signal-to-noise ratio for the network of the three independent interferometers is then 

1/2 

(26) 



P 



La=i 

In this paper, we shall focus on the estimation of the parameter lntf/,. The 1-a observational 
error, <r , can be estimated from the Fisher matrix Ay. An important point is that shGRBs are 
believed to be beamed: the 7 radiation is emitted in a narrow cone more or less perpendicular to 
the plane of the inspiral. We will take the total beaming angle to be at most 40° [32| (corresponding 
to t < 20°). It will be assumed that shGRBs are produced by the mergers of the neutron star 
binaries. For definiteness we take them to have component masses of (1.4, 1.4)M sun . We will 
consider 1000 sources up to a redshift of z = 2, which is where the angle-averaged signal-to-noise 
ratio approximately reaches the value 8 for sources with 1 < 20°. 

Before continuing, we mention that coalescing binaries composed of a neutron star and a black 
hole (NSBH) could also cause shGRBs 32|. For a fixed distance, a NSBH event will have a larger 



SNR than a BNS event, leading to an improved measurement of lnd^. The intrinsic event rates for 



NSBH are quite uncertain, but they are expected to be considerably lower than for BNS 33]; on 
the other hand, NSBH events will be visible to ET out to redshifts of z ~ 4 [231 ] . It is likely that the 
inclusion of NSBH would have a noticeable beneficial effect on the determination of cosmological 
parameters, especially if the black holes have spin, but this we leave for future studies. 

For a given event, distance measurements will be subject to two kinds of uncertainties: the 
instrumental error a a which can be estimated using a Fisher matrix as discussed above, and an 
additional error o\ due to the effects of weak lensing. As in previous work [23] we assume the con- 
tribution to the distance error from weak lensing to satisfy 07 = 0.05z. Thus, the total uncertainty 
on A In (if, is taken to be 



&hid L = ^o* + of. (27) 

In the next subsection we discuss how the information from multiple GW standard sirens can 
be combined to compute the expected measurement uncertainties on cosmological parameters. 

C. Gravitational- wave standard sirens 

Now let us turn to the determination of the cosmological parameters, including dark energy 
parameters, by the GW standard sirens. For each shGRB source, the luminosity distance is 
measured from the GW observation, and the redshift z can be obtained from the electromagnetic 



counterpart. Thus, the di — z relation can be employed to constrain various cosmological pa- 
rameters. For the cosmological model introduced in Sec. Ill Al we consider five free parameters 
(wo, w a , O m , fifc, ho) which can be constrained by GW standard sirens. We note that the value of 
tide (the relative energy density of the dark energy component) is determined by the quantities Cl m 
and through Q^e = 1 — Q m — ^k- 

In order to estimate the errors on these parameters, we study a Fisher matrix i ? i ^ w for a 
collection of inspiral events: 



^ (A\nd L (k)) 2 ' { ' 

k 

where the indices i and j run from 1 to 5, denoting the free parameters (wq, w a , Q m , fi^, ho). 
Eq. ()10p gives the expression for di, and the partial derivatives with respect to the parameters are 
evaluated at the parameter values corresponding to the fiducial cosmological model of Sec. Ill Al 
The uncertainty A lntf^ is calculated by using Eq. (|27p . The index k = 1, 2, . . ., labels the event at 
{ z kilk)i where the vector 7 stands for the angles (0, (ft, 1, tp). Here we should mention that, in (|28p . 
we have ignored the photometric redshift errors and the possible errors generated by the peculiar 
velocities of the sources relative to the Hubble flow [34j ]. Given that the majority of our sources 
will be at z > 0.4, we do not expect the latter to make much difference to our main results. 
Since d^ is independent of 7, this Fisher matrix can be rewritten as 

ft =x>cta<M*»%<ta«w%»{x; iSEazbsjp) • (29) 

z k \ 7fc ) 

When the number of events is large, the sum over events in (|29p can be replaced by an integral, so 
that we obtain 

F W = r O i (lnd L )d j (\nd L )f(z)A(z)dz, (30) 
J 

where f(z) is the number distribution of the GW sources over redshift z. A(z) is the average of 
1/(A In di{z, 7)) 2 over the angles (0, <j>, t, ip) with the constraint 1 < 20°: 

^H(Aln^(z,7)) 2 ), ; l<200 - (31) 
In order to calculate the averaged quantity A(z), we used a Monte Carlo sampling with 10,000 
choices of 7 for a given z, where z ranges from to 2 in steps of 0.1. The results are indicated in 
Fig. Q] by the red dots. We find that these points can be fit very accurately by a simple relation 
(see the black solid line in Fig. [T|) : 

A-V 2 (z) = 0.1449z - O.OII82 2 + 0.0012z 3 , (32) 
10 




FIG. 1: The averaged quantity A 1 ^ 2 {z) [defined in (|3Tj) ] as a function of redshift z. The red dots denote 
the results based on the Monte Carlo sampling, and the solid line denotes the fit results. 

which is used in our subsequent calculation. 

In (|30p . the upper integration limit z = 2 is the redshift at which the angle- averaged signal-to- 
noise ratio is approximately 8 [^J. In our fiducial model with $7/% = 0, the number distribution 



f(z) is given by 



_ 4nMr(z)d 2 c (z) 
!{Z) ~ H(z)(l + z) ' {66) 



where dc is the comoving distance, which is defined as dc{z) = Jq 1/ H(z')dz' . The function r(z) 
describes the time evolution of the burst rate, and the constant N (the number of the sources 
per comoving volume at redshift z = over the observation period) is fixed by requiring the total 
number of the sources -/Vgw = Jq f( z )dz. The expected total number of inspirals per year within 
the horizon of ET is ~ several x 10 5 for neutron star binaries. If, as suspected, neutron star binaries 
are progenitors of shGRBs 32], it might be possible to make a coincident detection of a significant 
subset of the events in the GW and electromagnetic windows, which can then be considered as 
standard sirens. As we have mentioned, shGRBs are believed to be beamed with small beaming 
angle, so only a small fraction of the total number of neutron star binaries are expected to be 
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1000 




FIG. 2: The normalized distribution of the GW sources 
observed as shGRBs. Following 



23l |. we assume that about 1000 events (~ 10 -3 of the total 
number of binary coalescences) will be observed in both windows, i.e., JVgw = 1000 throughout 
this paper. 

Since the time evolution of the source rate is as yet unclear, in this paper we shall consider two 
different forms for the function r{z). In the first case we assume that the sources are distributed 
uniformly, i.e., with constant comoving number density throughout the redshift range < z < 2 
(hereafter we will refer to this as the uniform distribution). In this case we have r(z) = 1, which 



is what was assumed in the previous work 



231 ] . In the other case, we take r{z) to be the following 



function: r{z) = (1 + 2z) for z < 1, r{z) = (15 — 3z)/4 for 1 < z < 5, and z = for z > 5. 



351 ] . Hereafter, we shall call this the 



This approximate fit to the rate evolution is suggested in 
nonuniform distribution. In Fig. [2l we plot the distribution function / as a function of redshift z 
in the two cases. Note that in the case with nonuniform distribution, the sources are a little more 
concentrated at z = 1. In what follows we will find out how this affects the uncertainties on the 
cosmological parameters. 

Using the definition (|30p . we can calculate the Fisher matrix i ? ^ w for the following two cases : 
a. GW events with uniform distribution; b. GW events with nonuniform distribution. The results 
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TABLE I: GW Fisher matrix in the case with uniform distribution 





W Q 


W a 


r2 m 


n k 


h 


w 


0.273681 x 10 4 


0.433681 x 10 3 


0.753068 x 10 4 


0.261216 x 10 4 


0.221885 x 10 5 




0.433681 x 10 3 


0.753710 x 10 2 


0.128638 x 10 4 


0.406152 x 10 3 


0.317462 x 10 4 




0.753068 x 10 4 


0.128638 x 10 4 


0.221212 x 10 5 


0.704345 x 10 4 


0.571447 x 10 5 


r» fe 


0.261216 x 10 4 


0.406152 x 10 3 


0.704345 x 10 4 


0.251558 x 10 4 


0.213013 x 10 5 


h 


0.221885 x 10 5 


0.317462 x 10 4 


0.571447 x 10 5 


0.213013 x 10 5 


0.216280 x 10 6 



TABLE II: GW Fisher matrix in the case with nonuniform distribution 





W Q 


W a 


r2 m 


n k 


h 




0.256794 x 10 4 


0.427648 x 10 3 


0.731269 x 10 4 


0.244368 x 10 4 


0.194634 x 10 5 


W a 


0.427648 x 10 3 


0.762633 x 10 2 


0.129200 x 10 4 


0.399934 x 10 3 


0.303753 x 10 4 




0.731269 x 10 4 


0.129200 x 10 4 


0.219941 x 10 5 


0.682599 x 10 4 


0.529628 x 10 5 




0.244368 x 10 4 


0.399934 x 10 3 


0.682599 x 10 4 


0.234666 x 10 4 


0.186267 x 10 5 


h 


0.194634 x 10 5 


0.303753 x 10 4 


0.529628 x 10 5 


0.186267 x 10 5 


0.162814 x 10 6 



are listed in Tables I and II. 

From the Fisher matrices, we can calculate the 1-a uncertainties on the parameters, which are 
A Pi = \J(F QW )~]i- For the case with uniform distribution, by using the results in Table I, we find 
that 

Au; = 2.62, Aw a = 9.53, Af} m = 0.815, AQ k = 2.03, Ah = 1.20 x 10~ 2 . (34) 

We plot the two-dimensional uncertainty contour of the parameters wo and w a in Fig. [3] (blue 
curve, i.e., line 2, in the left panel). Unfortunately, we find that the error bars on the parameters 
are all fairly large, especially for the dark energy parameters wq and w a . This is caused by the 
strong degeneracy between (wq, w a ) and the other parameters (O m , ho), especially (£l m , fife). 
To illustrate this, let us do the following calculation. First we fix the values of the parameters 
(fi m , fifc, ho) to be their fiducial values, and only consider (wo, w a ) as free parameters. By using 
the results in Table I, we obtain 

Aw = 0.064, Aw a = 0.388. (35) 

We find that the values of Awo and Aw a become much smaller in this case. The two-dimensional 
uncertainty contour of wo and w a is also plotted in Fig. [3] (black curve, i.e. line 1, in the left panel). 
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This figure shows that there is correlation between the parameters wq and w a . Recall that a goal of 
the dark energy programs is to test whether dark energy arises from a simple cosmological constant, 
(wq = —1, w a =0). For a given data set we can do better (as far as excluding the cosmological 
constant model is concerned) than simply quoting the values of Awo and Aw a . This is because the 
effect of dark energy is generally not best constrained at z = 0. For the phenomenological form of 
the EOS of the dark energy w(z) = wo+w a z/(l+z), the constraint on w(z) varies with the redshift 
z. So, similar to 6j, we can define the best pivot redshift, denoted as z p , where the uncertainty 
in w(z) equals the uncertainty in a model that assumes w a = 0. In this paper, we denote the 
EOS at this best pivot redshift as w p = w(z p ). The best pivot redshift z p can be calculated by 
z p = —1/(1 + f^ Q ), where p is the correlation coefficient of wq and w a . The value of Aw p is 
calculated by Aw p = Awq \J\ — p 2 . In this case (two free parameters) , the results for z p and Aw p 



arc 



z p = 0.188, Aw p = 0.019. (36) 

The value of Aw„ as well as that of Aw a are commonly used to describe the detection ability of 
the experiments |6|. 

On the other hand, we can also fix the values of the parameters (wq, w a ) to be their fiducial 
values, and only consider (O m , O^, ho) as free parameters. By using the results in Table I, we 
obtain 

An m = 0.021, AQ k = 0.087, Ah = 5.48 x 10~ 3 . (37) 

Again we find that the values of these errors, especially the values of AO m and Afifc, are much 
smaller that those in Eq. (|34p . These results show that the GW standard sirens can constrain the 
dark energy parameters rather well, on condition that we can break the strong degeneracy between 
the parameters (wq, w a ) and the parameters (fl m , fifc, ho). In the next subsection, we will find 
that this can be realized if we consider the CMB observations as a prior. 

To conclude this subsection we discuss the determination of the dark energy parameters (wo, 
w a ) by ET observations in two cases considered in this paper. By using the results in Table I and 
II, we calculate the errors of the two parameters (the other parameters are fixed at their fiducial 
values). The two-dimensional uncertainty contours are shown in Fig. [3] (upper left panel). This 
figure shows that the errors of the parameters are a little larger for the nonuniform distribution 
than those in the corresponding case with the uniform distribution. This is because the sources 
in the nonuniform distribution are a little more concentrated at the redshift z = 1. Hence the 
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FIG. 3: The two-dimensional uncertainty contours of the dark energy parameters wq and w a in the case 
with uniform distribution. 

number of the sources in the high redshift and the low redshift regions is smaller, making it more 
difficult to constrain the dark energy's evolution. So, in addition to the number of the sources, the 
redshift distribution of the sources also plays a crucial role for the detection of dark energy. The 
results listed in (|35l) constitute the optimistic case among the cases we have considered. On the 
other hand, it is helpful to list the results in the case with nonuniform distribution, which are 

Aw = 0.077, Aw a = 0.445. (38) 

We find that these uncertainties are a little larger than those in (1351) . We can also calculate the 
values of z p and Aw p for this case, which are 

z p = 0.200, Aw p = 0.020. (39) 

The value of Aw p is also a little larger than that in the case with uniform distribution. 

D. Planck CMB prior 

As will be clear from the discussion above, the ability of GW sources to constrain dark energy 
depends strongly on how well the background parameters Q m and Qk can be measured beforehand. 
Indeed, if one tries to determine these background parameters as well as the dark energy parameters 
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FIG. 4: The two-dimensional uncertainty contours for the dark energy parameters wo and w a , in the different 
cases. 

together using GW sources, the values of Au;o and Aw a become very large, and the constraints on 
dark energy become meaningless. Hence we should consider another observational method which 
can determine these background parameters through a prior observation. As we shall see, this is 
also necessary for the other ways used to study dark energy - BAO and SNIa - so that it does 
not diminish the value of having self-calibrating standard sirens: GW observations will provide us 
with an important independent check. 
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Observations of the Cosmic Microwave Background (CMB) temperature and polarization 
anisotropies are always used as the required prior. The WMAP satellite has already given fairly 
good results for the CMB TT (temperature-temperature auto-correlation) and TE (temperature- 
polarization cross correlation) power spectra in the multipole range I < 800. The background 
parameters (fl m , fife) have already been well determined by 7- year WMAP observations; for exam- 
ple, the constraint on the curvature is —0.0133 < Oj. < 0.0084 (95% C.L.) ^]. These constraints 
are expected to be greatly improved by the Planck observations [36(, which will give good data 
on the CMB TT and TE power spectra up to £ ~ 2000. In addition, Planck is also expected to 
observe the EE (electric type polarization) power spectrum. In this subsection, we shall consider 
the potential CMB observations by the Planck satellite. For the noise power spectra of Planck, we 
consider the instrumental noises at three frequency channels: at 100GHz, 143GHz, and 217GHz. 
For the reduced foreground radiations (including dust and synchrotron), we assume that the re- 
duced factor <r fg = 0.1, i.e. 10% residual foregrounds are considered as the noises. The total noise 
power spectra Nj T and N EE of the Planck satellite can be found in [37]. Note that in this paper 
we assume 4 sky (28 months) survey. 

In order to estimate the constraints on the cosmological parameters, we shall again use the 
Fisher information matrix technique. The Fisher matrix is calculated by 381 ] 

1=2 XX', YY' ° Pl ° P 3 

where pi are the cosmological parameters to be evaluated. C xx are the CMB power spectra 
and D xx are their estimates. Cov -1 is the inverse of the covariance matrix. The non-vanishing 
components of the covariance matrix are given by 

Cov(D xx ,D xx ) = {2i+ \ )fsky ( C ? X + N ^ 2 (X = T,E), 
Cov{Dj E ,Dj E ) = ^ + \ [{Cj E f + (Cj T + Nj T ){Cf E + AO], 

^(DF,Dr) = J ^^Cj E i C F + N F) , 
Cov(B? E ,Df E ) = {2e+ \ )fsky C ? E ( C f E + N z EE )- 
Note that in the calculation we have adopted ^ max = 2000, and the sky-cut factor f s ^ y = 0.65 



suggested by Planck Bluebook 



36j. 
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TABLE III: CMB Fisher matrix 





w 


W a 


r2 m 


n k 


ho 




0.414303 x 10 5 


0.115085 x 10 5 


0.287229 x 10 6 


-0.678690 x 10 6 


0.339923 x 10 6 


W a 


0.115085 x 10 5 


0.320204 x 10 4 


0.797415 x 10 5 


-0.190373 x 10 6 


0.943498 x 10 5 




0.287229 x 10 6 


0.797415 x 10 5 


0.219854 x 10 7 


-0.465663 x 10 7 


0.251813 x 10 7 




-0.678690 x 10 6 


-0.190373 x 10 6 


-0.465663 x 10 7 


0.136912 x 10 8 


-0.548091 x 10 7 


ho 


0.339923 x 10 6 


0.943498 x 10 5 


0.251813 x 10 7 


-0.548091 x 10 7 


0.291587 x 10 7 



The CMB power spectra C^ x ' depend on all the cosmological parameters, including the back- 
ground parameters and the perturbation parameters. In the calculation, we first build the Fisher 
matrix for the full nine parameters (wo,w a ,Q.i ) ti^,Q, c h^,Q<k,n s , A s ,ho,r). In order to obtain the 
constraint on the parameter O m , we change the full Fisher matrix to the new one of the nine 
parameters (wq, w a , ^ c ^0' ^k, n s, A s , ho, r), where Q m = (O&Ziq + f2 c ^o)/^o nas been used. In 
order to directly compare and combine with the Fisher matrix of GW method, we marginalize 
the new nine-parameter Fisher matrix to the one with five parameters (wq, w a , O m , fi^, ho). The 
results are shown in Table III. The errors of the parameters are given by Api = (F CMB )" 1 ii . By 
using the Fisher matrix in Table III, we obtain 

Aw = 0.411, Aw a = 0.517, Aft m = 8.88 x 10~ 2 , Att k = 2.27 x 1(T 3 , Ah = 0.115. (41) 

This result shows that Planck alone can give quite tight constraints on D, m and f2fc, which is just 
complementary with the GW method. However, CMB observations alone cannot constrain the 
parameters wq and w a , which is because that the CMB power spectra are only sensitive to the 
physics in the early Universe at z ~ 1100, where dark energy is totally subdominant. 

We now investigate the combination of CMB and GW methods. In order to do this, we define 
a new Fisher matrix, which is the sum of F^- and F^j MB . By using this new Fisher matrix, we 
obtain 

Aw = 0.079, Aw a = 0.261, Af2 m = 5.14 x 10 -3 , AfL k = 6.66 x 10~ 4 , Ah = 5.96 x 1(T 3 , (42) 

for the uniform distribution case. We find that the values of Awo and Aw a are fairly close to those 
in Eq. (|35p . where we have only considered the GW observations but assumed the backgrounds 
parameters Q m , fifc and ho are fixed. These results show that taking the CMB observation as a 
prior is nearly equivalent to fixing the parameters f2 m , Q k an d ho. The two-dimensional uncertainty 
contour of the parameters wo and w a is shown in Fig. [3] (red lines, i.e. line 3, in the left panel and 
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the middle panel). We can also calculate the best pivot redshift z p and the value of Aw p , which 
are 



By the same method, we also obtain the results in the nonuniform distribution case. The two- 
dimensional uncertainty contour of wo and w a is shown in Fig. [3] (upper-right panel). In this case, 
the results are 



Aw = 0.099, Aw a = 0.302, An m = 7.30 x 10~ 3 , AQ k = 6.70 x 10" 4 , Ah = 8.99 x 10" 3 . (44) 
The best pivot redshift z p and the uncertainty Aw p are 



Again as expected, we find that the values are larger than the corresponding values in the uniform 
distribution case. 

III. DETECTION OF DARK ENERGY BY BAO AND SNIA OBSERVATIONS, AND THE 
COMPARISON WITH ET-GW OBSERVATIONS 

In the above we found that by combining the potential Planck CMB observation with ET-GW 
observations, one can get fairly tight constraints on the dark energy parameters wo and w a . In this 
section we discuss the detection abilities of the other two probes: BAO and SNIa. Currently these 
two methods play the crucial role for the determination of the dark energy component. In the 
near future, the detection abilities of these two methods are expected to be significantly improved; 
in this section a detailed discussion is given. All three probes, BAO, SNIa and GW, constrain 
the EOS of dark energy by probing the large-scale background geometry of the Universe (different 
from the weak gravitational lensing method [39]), so a fair comparison can be made, as we shall 
do here. 

A. Detection of dark energy by potential BAO observations 

The BAO method relies on the distribution of baryonic matter to infer the redshift-distance 
relation. The characteristic scale length of structure which can be accurately determined from 
the CMB is used as a standard rod. By measuring the angular size of this characteristic scale- 
length as a function of redshift, the effect of dark energy can be inferred. The BAO method can 



z p = 0.401, Aw p = 0.025. 



(43) 



z p = 0.454, Aw p = 0.030. 



(45) 
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constrain the dark energy by two observable quantities In (H(z)) and In (cIa(z)), where H(z) is the 
binned Hubble parameter and cLa{z) is comoving angular diameter distance, which is related to 
the luminosity distance by cIa = cIl/(1 + z). Similar to the quantity di, these two observables only 
depend on the cosmological parameters wq, w a , £l m , &k and ho, which will be considered as the 
parameters determined by the observations. 

In order to investigate the constraints on the cosmological parameters, we build the following 
Fisher information matrix Qj: 

pBAO = y dln(H(z k ))dln(H(z k )) ( 1 \ 2 + 8 In (d A (z k )) 8 In (d A (z k )) ( 1 \ 2 (4g) 

13 ~ dPi dPj V°ln(H(*i))/ d Pi d Pj \0\n{d A { Zi )) ) 

Again the index k denotes the observables, which are binned into several redshift bins, pi denotes 
the cosmological parameters. (T\ n (H(z)) and &\n(d A (z)) are the errors (including the observational 
errors and the systematic errors) of the observables In {H{z)) and In (<1a(z)), respectively. In order 
to study the detection ability of the BAO method, we shall consider the potential observations of 
a typical project, the final JDEM (Joint Dark Energy Mission) project (J, which is expected to 
survey 10000 deg 2 in the redshift range z £ (0.5,2). In the calculation, we bin the observables 
\n(H{z)) and \n(dA(z)) into 10 redshift bins, i.e. Az = 0.15 for each bin. The calculation of the 
theoretical values of these quantities are straightforward. For the errors of these observable data, 
we use the fitting formulae derived in {^(J (see also Eq. (4.8) in [g]). 

The results of the Fisher information matrix are shown in Table IV. To begin with we consider 
a simple case with only two free parameters (wq, w a ). We assume that the other parameters (fi m , 
fifc, ho) are fixed to their fiducial values. By using the Fisher matrix in Table IV, we obtain the 
uncertainties of the free parameters: 

Aw = 0.087, Aw a = 0.346, z p = 0.323, Aw p = 0.023. (47) 

However, if we try to constrain all five parameters by BAO observations, the uncertainties of 
the parameters will become fairly large. For instance, the uncertainties of too and w a become 
Au>o = 0.850 and Aw a = 3.611, respectively, which are much larger than those in (|47p . Similarly 
to the discussion in Sec. Ill Dt we can consider the combination of the BAO observation and the 
Planck CMB prior. By analogous steps we obtain the results 

Aw = 0.176, Aw a = 0.415, An m = 2.01 x 10~ 2 , AQ k = 6.40 x 10 -4 , Ah = 2.57 x 10~ 2 . (48) 

The best pivot redshift and the value of Aw p are 

z p = 0.664, Aw p = 0.059. (49) 
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TABLE IV: BAO Fisher matrix 





W Q 


W a 


r2 m 


n k 


h 




0.193253 x 10 4 


0.471352 x 10 3 


0.865506 x 10 4 


0.288000 x 10 4 


0.125880 x 10 5 




0.471352 x 10 3 


0.123299 x 10 3 


0.234243 x 10 4 


0.788613 x 10 3 


0.309758 x 10 4 




0.865506 x 10 4 


0.234243 x 10 4 


0.457760 x 10 5 


0.153244 x 10 5 


0.579934 x 10 5 


r» fe 


0.288000 x 10 4 


0.788613 x 10 3 


0.153244 x 10 5 


0.541597 x 10 4 


0.190452 x 10 5 


h 


0.125880 x 10 5 


0.309758 x 10 4 


0.579934 x 10 5 


0.190452 x 10 5 


0.834616 x 10 5 



B. Detection of dark energy by potential SNIa observations 

Now, let us turn to discuss the detection of dark energy parameters by the SNIa probe. Type 
la supernovae serve as a standard candle of (approximately) known luminosity. The redshift of 
supernova can be obtained by studying its spectral lines. Thus the redshift-distance relation can 
be gotten from SNIa surveys. Now, SNIa observed from various experiments have been used 
successfully to deduce the acceleration of the Universe after z = 1 jj, |4l(. In the near future, 
observations of SNIa are expected to be significantly improved, so that they will continue to serve 
as one of the most important methods for the determination of dark energy. 

The observables for SNIa data are the apparent magnitudes m, which can be corrected to behave 
as standard candles with absolute magnitude M with m = M + /j,(z). The function /i(z) for the 
measured redshift is 

fj,(z) = 5log 1Q (d L (z)) + 25, (50) 
where di is the luminosity distance. In this paper, we shall consider SNIa observations by the 



future SNAP (Supernova/Acceleration Probe) project 421 ] . As suggested by the SNAP group, we 
consider 300 low redshift supernovae, uniformly distributed over z £ (0.03,0.08). The error bar 
on the magnitude is assumed to be a m = 0.15 mag. In addition, 2000 high redshift supernovae in 
the range z £ (0.1, 1.7) are considered. The expected redshift distribution of these sources can be 
found in the SNAP white book (the middle red curve in Fig. 9 of [42| ) . We bin these 2000 sources 
into 10 redshift bins in the range z £ (0.1, 1.7). The total errors of the observables &{z) can be 
estimated as follows |3] : 



a = ^al + al, (51) 
where a\ = 0.15 mag/V^ (N is the total number of supernovae in each bin) is the intrinsic random 

21 



TABLE V: SNIa Fisher matrix 





W Q 


W a 


r2 m 


n k 


h 


w 


0.703955 x 10 4 


0.122773 x 10 4 


0.207309 x 10 5 


0.669019 x 10 4 


0.522067 x 10 5 




0.122773 x 10 4 


0.225085 x 10 3 


0.377761 x 10 4 


0.115185 x 10 4 


0.849589 x 10 4 




0.207309 x 10 5 


0.377761 x 10 4 


0.636028 x 10 5 


0.194262 x 10 5 


0.146902 x 10 6 


r» fe 


0.669019 x 10 4 


0.115185 x 10 4 


0.194262 x 10 5 


0.639862 x 10 4 


0.498026 x 10 5 


h 


0.522067 x 10 5 


0.849589 x 10 4 


0.146902 x 10 6 


0.498026 x 10 5 


0.491151 x 10 6 



Gaussian error, and 02 = 0.02mag(l + z)/2.7 is the error due to the astrophysical systematics. 
Thus, we can build a Fisher information matrix, which is 

f sn _ dfi(z k ) dfi(z k ) ( 1 \ 2 

The results are shown in Table V. The errors of the cosmological parameters are estimated by 
Api = J (-F SN ) _1 jj- Similar to Sec. IIII Al we first discuss the simplest case with two free parameters 
(wo, w a ). The results are 

Aw = 0.054, Aw a = 0.302, z p = 0.211, Aw p = 0.012. (53) 

To constrain all five cosmological parameters, we consider the combination of SNIa and the Planck 
CMB prior to decouple the degeneracy between (wq, w a ) and the other parameters. The errors of 
the parameters are 

Aw = 0.051, Aw a = 0.201, AQ, m = 3.49 x 10" 3 , AVt k = 6.52 x 10~ 4 , Ah = 3.39 x 10~ 3 . (54) 

We find that the values of Awo and Aw a in ([54|) are close to those in (i53|) . The best pivot redshift 
and the value of Aw p are also obtained 

z p = 0.313, Aw p = 0.019. (55) 



C. Comparison with the ET-GW observations 

Now let us compare the detection abilities of various probes: GW, BAO and SNIa. First we 
consider the simplest case, where only dark energy parameters (wq, w a ) are considered. The errors 
of the parameters are given in (I47p for BAO, and in (153p for SNIa. We find that the values of Au>o, 
Aw a and Aw p are all smaller for the SNIa probe. This shows that, comparing with the JDEM 
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BAO project, the SNAP SNIa project is expected to give a tighter constraint on the dark energy. 
For the ET-GW project, in Sec. Ill C( we have considered two cases. For the uniform distribution 
case, the results are given in (|35|) and (|36p . On the other hand, if the nonuniform distribution case 
is considered, the results are shown in ([38]) and ([39]) . Comparing with these results, we conclude 
that, in both cases, the detection ability of ET-GW method is stronger than that of the JDEM 
BAO project, but weaker than that of the SNAP SNIa project. This is mainly because the number 
of the GW standard sirens (~ 1000 as we have assumed) is smaller than that of the SNIa standard 
candles (~ 2000 at high redshift). In order to clearly show this, let us consider another case for ET- 
GW method, where we assume that 2000 sources with nonuniform distribution will be observed. 
We then obtain Awo = 0.054, Aw a = 0.315 and Aw p = 0.014, which comes close to the projected 
uncertainties of the SNAP SNIa project given by (|53|) . We note that the relative error bars of the 
SNIa in Eq. (51) are larger than those of GW sources in (32), especially at low redshifts. However, 
we find that this disadvantage of the SNIa method is overcome by the assumed 300 low redshift 
SNIa in z £ (0.03,0.08). 

We can also compare the results of these three probes, when considering the full 5 cos- 
mological parameters and adopting the Planck CMB prior. In Fig. [3] (middle panel), we 
plot the two-dimensional uncertainty contours of the parameters (wq, w a ), where for the ET- 
GW method we have considered the case with uniform distribution. Similarly, we find that 
the red curve (GW+CMB) is only a little looser than the green one (SNIa+CMB), but 
much tighter than the black one (BAO+CMB). In Fig. [3] (right panel), we plot the re- 
sults of the two-dimensional uncertainty contours for the four combinations (CMB+BAO+SNIa, 
CMB+GW+BAO, CMB+GW+SNIa, GW+SNIa+BAO). We find that the first three combina- 
tions have similar results. However, the constraint on the dark energy parameters is much looser 
for combination of GW+SNIa+BAO, where the Planck CMB probe is absent. This panel shows 
that the CMB prior indeed plays a crucial role in the detection of dark energy. 

Let us now combine all the four probes to constrain the cosmological parameters, including the 
dark energy parameters. If we consider the uniform distribution for the GW sources, we obtain 
the constraints on the dark energy parameters 

Aw = 0.044, Aw a = 0.171, z p = 0.308, Aw p = 0.017. (56) 

This is the best constraint what we could expect to obtain. If we consider the case with nonuniform 
distribution for the GW sources, the constraints slightly loosen to 

Aw = 0.045, Aw a = 0.174, z p = 0.313, Aw p = 0.017. (57) 
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In Fig. U] (lower-right panel) , we plot the two-dimensional uncertainty contours of the parameters 
(wq, w a ) for both cases. We find that the two curves are very close to each other; the relative 
weight of the GW probe is not very high for the combined methods. 

Finally, we would like to know how much the ET-GW probe can contribute in constraining all 
the five cosmological parameters (wq, w a , m , f^, ho). In order to do so, we first calculate the 
constraints on the parameters by the combination of CMB+BAO+SNIa. We obtain 

Aw = 0.048, Aw a = 0.184, AVt m = 3.46 x 10~ 3 , AQ k = 5.91 x 10~ 4 , Ah = 3.36 x 10~ 3 . (58) 

If we then add the contribution of the ET-GW probe with nonuniform distribution (the uniform 
distribution case gives the very close results), the results become 

Aw = 0.045, Aw a = 0.174, An m = 3.39 x 10" 3 , AQ k = 5.83 x 10" 4 , Ah = 3.20 x 10~ 3 . (59) 

In this case we find that, due to the contribution of ET-GW probe, Awo is decreased by a 6.3% and 
Aw a by 5.5%. In Fig. [5] we plot the two-dimensional uncertainty contours of the free parameters 
(wq, w a , Q m , fifc, ho) for both combinations. 



IV. CONCLUSION 



If short-hard 7-ray bursts (shGRBs) are produced by the mergers of neutron star binaries, the 
luminosity distances di of the sources can be determined by the Einstein Telescope gravitational- 
wave detector in the redshift range z < 2. The redshifts z of the sources can be determined with 
great accuracy through their electromagnetic counterparts. Thus it will be possible to use shGRBs 
as 'standard sirens' to study the dark energy component in the Universe by determining the EOS 
and its time evolution. 

When calculating the uncertainties in the determination of di by the ET observations, we 
assumed that the 7-ray emission is confined to a cone with an opening angle as large as 40° 



321 ]. corresponding to inclination angles 1 < 20°. In order to study the effect of the redshift 
distribution of the sources, we considered two different kinds of distributions: one in which sources 
are distributed uniformly in comoving volume, and a nonuniform distribution as in 3^]. We found 
that, by taking into account the Planck CMB prior, the errors on the dark energy parameters are 
expected to be Au>o = 0.079 and Aw a = 0.261 in the uniform distribution case, which is close 
to the detection ability of the SNAP Type la Supernovae project. Even in the 'pessimistic' case 
with nonuniform distribution of the sources, the errors are Awq = 0.099 and Aw a = 0.302, which 



24 




FIG. 5: The two-dimensional uncertainty contours of the cosmological parameters. The black (larger) 
curves denote the results of CMB+BAO+SNIa, and the red (smaller) curves denote the results of 
CMB+BAO+SNIa+GW. Here, we have assumed the GW with nonuniform distribution (The case with 
uniform distribution gives the nearly overlapped results ). 



is weaker than the detection ability of the SNAP Type la Supernova project, but stronger than 
that of the JDEM Baryon Acoustic Oscillation project. We also found that, comparing with the 
combination of the future CMB(Planck)+BAO(JDEM)+SNIa(SNAP) projects, the contribution of 
this kind of standard sirens can decrease the error of u>o by ~ 6.3% and that of w a by ~ 5.5%. Thus, 
the kind of self-calibrating GW standard sirens accessible to Einstein Telescope would provide an 
excellent probe of the dark energy component. 

Finally, it is important to mention that, in addition to GW, CMB, BAO and SNIa methods, 
there are a number of other probes, including cosmic weak lensing, galaxy clustering, and so 



on, which can also be used to detect the dark energy component in the Universe (see [6j, [44] 



for details). In practice, one should combine all the probes. However, in this paper, we have 
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emphasized that using shGRBs as 'standard sirens' constitutes an important complement to the 
general electromagnetic methods. 
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